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We estimate the dipole radiation emitted in fusion pro- 
cesses. We show that a classical bremsstrahlung approach 
can account for both the preequilibrium and the thermal 
photon emission. We give an absolute evaluation of the pre- 
equilibrium component due to the charge asymmetry in the 
entrance channel and we study the energy and mass depen- 
dence in order to optimize the observation. This dynamical 
dipole radiation could be a relevant cooling mechanism in the 
fusion path. We stress the interest in experiments with the 
new available radioactive beams. 

PACS numbers: 24.30.Cz; 25.70.Jj; 25.70.Lm; 25.60.Pj 
In the early stage of a heavy ion collision large ampli- 
tude collective motions can be excited. Monopole, dipole 
or quadrupole modes are typical examples. Tracing in- 
formation from their decay provides the opportunity to 
learn about the features of fermionic systems in particu- 
lar, extreme conditions. 

If the colliding nuclei have different N/Z ratios the 
charge equilibration process takes place. Several works 
have suggested that the related neutron-proton collec- 
tive motion has the features of a Giant Dipole Resonance 
(GDR) mode [00. In a microscopic approach based on 
semiclassical BNV transport equations ||, it has been 
studied in detail how a collective GDR responsenbde- 
velops in the entrance channel. In connection to this 
analysis three main phases can be identified. First, dur- 
ing the approaching phase the two partners, overcom- 
ing the Coulomb barrier, still keep their own response. 
Then it follows a dinuclear phase when the conversion of 
relative motion energy in thermal motion starts to take 
place, mainly due to nucleon exchange. The composite 
system is not thermally equilibrated and manifests, as a 
whole, a large amplitude dipole collective motion. The 
charge equilibration has initiated. Finally, during the 
third phase, the Compound Nucleus ( CN) formation, the 
nucleus is thermally equilibrated. In a phonon approach 
P,|To[ it was shown || for the first time that an enhance- 
ment of the GDR gamma emission should be observed 
if the number of GDR phonons, n^^^ when the CN is 
formed, is larger than the value corresponding to a sta- 
tistical equilibrium of the GDR (as oscillator) with the 
CN (as heat bath). This effect has been experimentally 
evidenced in several fusion reactions pd|-p^|. 

However it is important to notice that also in the 
second (dinuclear) phase a GDR is excited on top of 
nonequilibrium states |Tj| (Brink- Axel hypothesis). The 
number of GDR phonons is even larger than tIqj- )R |^] 
and consequently some pre-equilibrium GDR photons 



can be emitted also at this stage. In the attempt to es- 
timate this contribution we cannot use straightforwardly 
the "phonon" approach |5|,^| which is based on the as- 
sumption of the existence of a thermally equilibrated nu- 
cleus. A way to face this problem is to directly apply 
a bremsstrahlung ("bremss" ) approach. Considering the 
evolution of the collective dipole acceleration from the 
time when it suddenly rises, at the beginning of the sec- 
ond phase, until it is completely damped to a pure "ther- 
mal" component, we can consistently calculate the whole 
contribution of the pre-equilibrium GDR to the photon 
yield. This is the aim of our letter. 

The total photon emission probability from the dipole 
mode oscillations, as given by the bremsstrahlung for- 
mula, can be expressed as jl(| (E 7 = hui): 
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where X"(u>) is the Fourier transform of the acceleration 
X"{t) associated with the distance between the centers 
of mass of protons (P) and neutrons (N), X = R p — R n . 
A = N + Z is the system mass. Thus following the time 
evolution of the dipole mode along the fusion dynamics 
it is possible to evaluate, in absolute values, the corre- 
sponding pre-equilibrium photon emission. However it is 
interesting to notice that the "bremss" approach [[l7) ac- 
tually provides a unified picture to account for the photon 
emission from both pre-equilibrium and thermal GDR. 
In fact Eq.(l) also allows to recover the usual expres- 
sion for the 7-decay rate from a GDR in thermal equi- 
librium with a CN at temperature T. In this case the 
dipole mode will manifest thermal oscillations that by 
the fluctuation-dissipation theorem are completely deter- 
mined only by its dissipative properties Jl8| . We may 
consider equivalently that these are the consequence of 
the action of a fictious, time dependent, random force 
F(t) |19|. Then the Fourier transform of the distance 
X(t) can be expressed as: X(u) = a(uj) x F(u>), where 
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is the response function of the dipole mode [j21| . M co u = 
^£-rn with m = 935MeV (nucleon mass) represents the 
collective mass of the neutron-proton relative motion and 
r is the decay width. 

The spectral density of mean square fluctuation of ran- 
dom force < \F(oj)\ 2 > is related to the properties of the 
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system through the fluctuation-dissipation theorem, as it 

follows raj 
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Using the Parseval theorem and the ergodic hypoth- 
esis we have, for times T very large, limj^oo T < 
\F(lu)\ 2 >= \F(uj)\ 2 . So, by introducing in Eq.(§) 
the Fourier transform X"(u>) of the acceleration, since 
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|X(w)|, we get the average photon emis- 



sion probability per unit time: 
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If we exclude the zero-point motion contribution and ex- 
pand the hyperbolic function for Klo > T we obtain: 
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where a a b s is the 7-absorption cross section in the GDR 
region. This is the well known formula for the gamma 
emission rate from hot GDR (^o|. It was previously ob- 
tained by using the detailed balance principle or from 
considerations on the statistical equilibrium between 
dipole mode and black-body radiation (2l| . 

Now we focus on the dynamical dipole mode in the 
entrance channel and we study the dependence of the 
pre-equilibrium dipole oscillations on the incident beam 
energy. We look at the relative weight vs. the statistical 
contribution. 

In the following we present results concerning the 
charge equilibration process for the systems 40 Ca 
(N/Z = 1) on 100 Mo (N/Z = 1.38), at AAMeV and 
w O (N/Z = 1) on 98 Mo (N/Z = 1.33) at various beam 
energies (4, 8, 14 and 20AM eV). Our fusion dynamics 
is obtained in the framework of an accurate BNV trans- 
port simulation of central collisions, see || for a detailed 
analysis. 

The time evolution of the dipole moment D(t) = 
~X{t) and of quantity DK(t) = TL/% for the Ca + Mo 
reaction is shown in Fig. la. Here II = — jf) with 

Pp (Pn) center of mass in momentum space for protons 
(neutrons) is just the canonically conjugate momentum 
of the X coordinate, II = M co uX'(t). Then a collec- 
tive GDR hamiltonian can be introduced as Hgdr = 
Ii 2 /(2M cM ) + M cM uj 2 X 2 /2. 

We choose the origin of time at the beginning of the 
dinuclear phase S . The dipole acceleration as obtained 
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by a second order numerical derivative of D(t) is shown 
in Fig. lb. With a good approximation this is zero be- 
fore t = 0. It suddenly rises when the collective dipole 
excitation is triggered and after few oscillations it be- 
comes completely damped, around 250/m/c. The corre- 
sponding power spectrum, \D"(oj)\ 2 , calculated from the 
Fourier transform of D"(t): 
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is plotted in Fig.lc by a dashed line. 
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FIG. 1. System Ca+Mo at 4 AMeV: (a) Time evolu- 
tion of D(t) (solid line, in fm units), and DK(t) (dashed 
line, in /m" 1 ). (b) The same for the acceleration D"(t) 
(in c 2 //m). (c) Power spectrum |D"(u>)| 2 (in c 2 units), (d) 
Bremsstrahlung spectrum (solid line) and the first step sta- 
tistical spectrum (dashed line). 

In the same Fig.lc we have added (solid line) \D"(lu)\ 2 
as obtained directly from the Fourier transform of the 
dipole moment, D(u>) by using: 

\D"(lo)\ 2 = D'(0) 2 + uj 2 D(0) 2 + 2uj 2 D'(0)ReD(uj) 

- 2uj z D(0)ImD(uo) + iv 4 \D(u;)\ 2 (7) 

where D(0) and D'(Q) are the dipole and the dipole ve- 
locity values at t — 0. The nice agreement between the 
two procedures indicates the numerical accuracy of the 
method. 

The total photon emission probability from the pre- 
equilibrium dipole mode, given by the bremsstrahlung 
formula Eq. (Q) , is reported in Fig. Id (solid line) . 

For comparison we show the first step statistical spec- 
trum (dashed line). The latter is just the product of 
the statistical gamma decay rate p0| , (see also Eq.(|h), 
times the mean life time of the compound nucleus tcn- 
A good estimate of tq n can been derived from the total 
neutron width of a CN at temperature T: T n — Ti/tcn = 
2mr 2 A 2 / 3 /(ttK 2 )T 2 exp(-B n /T), with B n = 8.5MeV 
(neutron binding energy) and ro = 1.2/m. For the sta- 
tistical GDR, the standard parameters from systemat- 
ics, including temperature dependence of the width, are 
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considered. We remark that the two contributions are 
comparable. 

The pre-equilibrium spectrum is shifted toward lower 
values of energy as a consequence of the large quadrupole 
deformation along the fusion path. However part of the 
total shift is also due to the used effective interaction M . 

We perform an analogous analysis for the reaction 1 °0 
+ 98 Mo at four beam energies. In Fig. 2 the power spec- 
tra \D"{tu)\ 2 obtained with the two procedures described 
before are plotted. Again the two methods give close 
results. 
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FIG. 2. Power spectrum \D"(uj)\ 2 for the system O+Mo 
at the beam energies 4 (a), 8 (b), 14 (c) and 20 (d) AMeV 
as obtained directly from dipole acceleration (solid line) and 
from Eq.(Q) (dashed line). 

In Fig. 3 the corresponding "bremss" spectra (solid 
lines) as well as the first step statistical spectra (dashed 
line) for the system O + Mo are plotted. We remark 
the characteristic "rise and fall" behaviour of the pre- 
equilibrium contribution (already observed within the 
"phonon" approach Q). This effect is more clearly ev- 
idenced in Fig. 4 where the extra photon multiplicities 
due to the charge asymmetry in the entrance channel, 
obtained integrating over energy in the resonance region, 
are plotted as a function of excitation energy for both 
"bremss" and "phonon" (see following) models. The re- 
sulting trend is determined by the interplay of several ef- 
fects. At low beam energies, a slow neck dynamics does 
not sustain the dipole oscillations. A strong attenuation 
of the acceleration takes place along entrance channel 
and the dipole emission is quite low. We expect a similar 
behaviour in deep-inelastic reactions. 

At higher beam energies a faster fusion dynamics can 
favor the dipole oscillation. Larger amplitudes of the ac- 
celeration and more dipole oscillations enhance the pre- 
equilibrium radiation. Around 8-14 AMeV an optimum 
effect is attained. At higher beam energies a larger GDR 
damping manifests and consequently the gamma emis- 
sion is again reduced || . 

From a comparison between the two systems, at the 
same energy available in the cm., we observe a larger 



pre-equilibrium emission from Ca + Mo (black square in 
Fig. 4). This can be related to larger values of the dipole 
acceleration as a result of a larger initial dipole moment 
amplitude. 

It is interesting to compare our results to the pre- 
dictions of the phonon approach, where the main in- 
put quantity is the number of GDR phonons at the 
time of CN formation. In this model a GDR phonon 
gas is coupled to the CN ||. The phonons are de- 
caying with a rate fi = T/h and excited with a rate 
A. The enhancement of 7-decay probability is given by 

APphonon = (nQDR ~ Vm^/G" + lev) were 7 7 is the 
partial width for photon emission and j ev is the evapora- 
tion decay rate. In this method one of the uncertainties 
is related to the estimation of the quantity nQ DR which 
in turn requires the knowledge of the CN formation time 
H . The results presented here are obtained following the 
usual procedures in r— and p— space to fix the thermal 
equilibration time, see ref. ||. We see from Fig. 4 that 
the values obtained in the "phonon" approach (circles) 
are systematically below the "bremss" evaluation, as ex- 
pected from the previous discussion. In order to get the 
"bremss" results we need much larger values of the pa- 
rameter n^p DR , corresponding to times well before statis- 
tical equilibrium, not consistent with the phonon model 
picture. 
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FIG. 3. Photon spectra for the O+Mo system at the four 
energies as in Fig. 2: bremsstrahlung (solid line) and first 
step statistical contribution (dashed line). 

In Fig. 4 we report also the ratio between total pre- 
equilibrium and total statistical photon multiplicity as 
a function of the initial excitation energy of the com- 
poud system. This ratio is a good measure of the pre- 
equilibrium effect which should be seen in experiments 
when the total photon yields from a N/Z asymmetric re- 
action and a symmetric one, both forming the same CN 
in similar conditions (temperature, angular momentum) 
are compared. To calculate the total statistical GDR 
contribution we consider a very schematic cascade Mon- 
tecarlo model including only gamma and neutron chan- 
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nel competition. The optimum range of beam energies to 
observe the pre-equilibrium effect appears to be around 
8AMeV. At higher beam energies a larger statistical 
contribution to gamma decay will lower the relative im- 
portance of pre-equilibrium contribution, which actually, 
as we have already discussed, starts to decrease in abso- 
lute value. We have to remark that, at the same E* , the 
effect is still more important for the Ca + Mo system, 
although now the total statistical emission is also larger, 
see Eq.(|^). The ~ 16% enhancement is also in very good 
agreement with the experimental estimate of pM- 



ics is quite fast. Exotic (radioactive) beams are a good 
opportunity for studying these phenomena. However a 
caution is related to the fact that now a charge equili- 
bration can take place also by an enhancement of fast 
nucleon emission. 

We finally remark that the dynamical nature of the 
pre-equilibrium contribution will show up in a clear 
anisotropic 7-emission since it is due to dipole oscillations 
on the reaction plane (see details in [pi). This could be a 
nice characteristic signature of fusion paths in reactions 
induced by radioactive beams. 
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FIG. 4. (Left) Extra photon multiplicity as a function 

of initial excitation energy: bremsstrahlung (empty squares), 

phonon model (empty circles) for the O + Mo system and 

Ca+Mo (the corresponding full symbols). (Right) Ratio of 

total pre-equilibrium to total statistical GDR 7-ray emission 

(see text). 

In conclusion, in this letter we have presented a 
bremsstrahlung approach to dipole radiation from both 
thermal GDR and dynamical dipole mode in fusion re- 
actions. We stress that in our calculation we directly 
get absolute values for the emission probabilities with- 
out any normalization or adjusted parameters. Moreover 
the "bremss" approach can be extended in a similar way 
to dissipative deep-inelastic collisions [0 where similar 
prompt 7 — GDR emission has been recently observed 

We have shown that when we consider the full dynam- 
ics of the pre-equilibrium dipole in charge asymmetric fu- 
sion processes the extra radiation contribution appears to 
be systematically two times larger than the expectations, 
based on a phonon model, used so far P ^j^l jTllJl^l (see 
Fig. 4). The pre-equilibrium dipole radiation due to the 
charge asymmetry in the entrance channel can be then 
comparable to the statistical CN emission. In this sense 
it can represent a new "cooling" mechanism in warm fu- 
sion reactions to form less excited residues, e.g. of inter- 
est for the synthesis of heavy elements. 

A dependence of the pre-equilibrium gamma emission 
on beam energy has been evidenced. We conclude that 
this is a consequence of the interplay between fusion dy- 
namics, damping properties of the dipole and statisti- 
cal gamma emission. We expect a better observation of 
the effect in reactions between nuclei with a larger initial 
dipole amplitude, at beam energies where fusion dynam- 
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